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ABSTRACT 



An algorithm has been designed to provide near real-time location 
and classification information to an operator of an underwater acoustic 
imaging system. The overall system to detect objects buried up to 
five meters in unconsolidated marine sediments consists of a spe- 
cialized projector co-located with a five meter receiving line array. 
This system is described in basic terms together with the principles of 
the signal processing that extracts information from the scattered 
acoustic signals. Using simulated objects, the angular resolution and 
imaging performance of a number of system designs were evaluated. To 
obtain near real-time execution, the algorithm was optimized and parti- 
tioned using parallel, pipeline, and double buffering techniques for 
independent but synchronized operation on multiple microprocessors. 

The feasibility of the design approach was experimentally demonstrated 
using four single board computers in a microcomputer development system. 
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I. INTRODUCTION 



The specific situation which led to this research is the desire 
to detect, locate, and classify objects buried in unconsolidated 
marine sediments up to a depth of 5 meters. To accomplish this an 
acoustic imaging system is needed which can be mounted on an under- 
water vehicle maneuvering close to the ocean floor. The problem is to 
penetrate the soft sediment acoustically, detect the scattered signals, 
and derive from them information as to the nature of the scattering 
objects. With sufficient range and angular resolution it would be 
possible to adequately localize and classify an object. Furthermore, 
it is necessary to obtain this information as rapidly as possible: 
ideally a display would provide a "real-time" representation of buried 
objects to an operator as a section of sediment is being searched 
acoustically. 

Although a number of techniques of acoustic imaging exist [.lj, 
a system similar to a high resolution sonar has been proposed which 
would have high resolution using very short baseband pulses and which 
would operate with dynamic focussing deep in the near field of the 
receiving aperture. This system would involve a specialized para- 
metric projector co-located with a receiving horizontal line array of 
hydrophones. It would be mounted on an underwater vehicle operating at 
about 5 meters above the ocean floor, and would provide acoustic 
penetration of the marine sediments to depths of about 5 meters. The 
objective of this thesis was to design the signal processing algorithm 
for such a system. 
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In the course of the development of the algorithm the effect of a 
number of system parameters on angular resolution and imaging perform- 
ance was evaluated for a number of possible implementations of the 
acoustic system. This led to major decisions as to the form of the 
algorithm and the possible acoustic imaging system. 

Since an acoustic imaging system mounted on a maneuvering under- 
water vehicle would be required to provide an operator with a near 
real-time display of buried objects, the algorithm was redesigned for 
increased execution efficiency and implementation on multiple micro- 
processors. Optimization of the algorithm by removing time consuming 
operations, together with incorporation of parallel and pipeline 
processes and double buffering concepts were the approaches taken. 
Experiments were conducted using a number of single board computers to 
test and demonstrate the feasibility of near real-time execution of the 
algorithm. 

Section II of the thesis describes the acoustic imaging system and 
the basis for the signal processing algorithm. The details of the 
algorithm and the angular resolution and imaging performance evaluations 
are explained in Section III. The design of the algorithm for near 
real-time execution on multiple microprocessors and the experimental 
demonstration of feasibility are contained in Section IV. The resulting 
conclusions are then given in Section V, followed in Section VI by 
recommendations for future work that could produce a special purpose 
processor for the acoustic imaging system. 
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II. THE PROPOSED ACOUSTIC IMAGING SYSTEM 



A. INTRODUCTION 

An acoustical Imaging system mounted on an underwater vehicle will 
consist essentially of three major parts: a sophisticated projector, a 

receiving array, and a signal processor and display system. The thrust 
of this thesis has been to develop the signal processing algorithm. To 
understand the underlying concepts, the potential form of the projector 
and receiving array is briefly described, followed by the proposed 
algorithm. 

B. TRANSMISSION 

Acoustic waves can penetrate sediment to useful depths with 
acceptable absorption loss if the frequency is low enough; less than 
about 20 kHz. Very short acoustic pulses (50 to 100 microseconds) are 
required for adequate range and angular resolution. For the concept 
that has been developed in this thesis it has become apparent that it 
would be necessary to have effectively a point source for the projector. 
Baseband pulse generation using a parametric source [2] appears to be a 
promising technique for this application. Such a projector would form 
a pencil beam in the collimated region of its near field that could be 
scanned to produce a tomographic image of a region in the sediment. 

C. RECEPTION 

1 . Physical Configuration 

The scattered acoustical energy would, be received by a horizon- 
tal line array (mounted on the underwater vehicle) of a large number of 
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equally spaced hydrophones. The projector would be located at the center 
of the line array. The length of the array would be a compromise between 
the need for a large aperture (on the order of 10 meters for angular 
resolution) and the physical reality of deploying and maneuvering an 
underwater vehicle with the array attached (5 meters probably the 
maximum). 

Moving over the area of interest the underwater vehicle would 
maintain an altitude of approximately 5 meters above the sediment. As 
the desired depth of penetration into the sediment is also 5 meters the 
total range of interest lies between 5 and 10 meters from the center of 
the array. A 60 degree sector (between array bearings of -30 and +30 
degrees from the normal to the line array) provides horizontal 
coverage of 10 meters in the sediment at 10 meters range from the 
center of the line array (Figure 1). 

2. Wavefront Curvature 

For ranges between 5 and 10 meters, operation of the acoustical 
imaging system is deep in the near field of the line array. The far 
field begins at a range of D^/A where D is the diameter of the aperture 
and A is the wavelength [1] . Using a 50 microsecond pulse and a 5 
meter line array this is on the order of 300 meters. The near field 
operation is evident by the marked curvature of the wavefronts of the 
acoustic energy as they arrive at the line array after scattering from 
points within the area of interest. Since the curvature depends on the 
range to a scattering point an imaging processor requires some form of 
dynamic focussing. 
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Figure 1. The physical relationship between the 
five meter line array of hydrophones 
and the sector of interest. 
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Figure 2 depicts four such wavefronts (that originated at the 
indicated bearings and ranges) as a function of time of arrival and 
position along a 5 meter line array. Using 1500 meters per second as 
the soundspeed, the time of the earliest arrival of interest is 6.220 
milliseconds after pulse transmission, and the latest is 14.304 milli- 
seconds . 

3. The “Trace 11 of a Point Scatterer 

The line array consists of a large number of equally spaced 
discrete hydrophones (each considered as a point) which detect the 
pressure amplitude of an arriving wavefront. The time sampling of the 
outputs of 16 hydrophones in a 5 meter line array resulted in a quanti- 
zation of the wavefront arrival times as shown in Figure 3 (in which 
an acoustic pulse width of 50 microseconds and sampling every 50 micro- 
seconds, beginning 6 milliseconds after pulse transmission, was used). 
For example, considering the wavefront originating at a scattering 
point located at a bearing of zero degrees and a range of five meters, 
the hydrophone sampling produces the following time indices for the 
wavefront arrivals: 21, 19, 17, 16, 15, 14, 13, 13, 13, 13, 14, 15, 

16, 17, 19, 21. 

The sequence of time indices representing the arrival of a 
wavefront, n_( r , 9 ) 3 (n^ >t\ 2 * . .n^) for K sensors, is unique for the 
scattering element at a range r and bearing 9: this sequence is termed 
the "trace" of that particular point (related to the concept of the 
"Time Delay Trace Function" [3] ). The traces for all points at the 
range and bearing increments in the area of interest can be pre-calcula- 
ted and stored for later retrieval . 
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Figure 2. The arrival times of four wavefronts as 

a function of position along a five meter 
array. 
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D. SIGNAL PROCESSING 



1 . The Basic Concept 

In developing the algorithm two situations were examined: 

a. Insonifi cation of the entire sector by a short pulse 
(broad beam in the plane of the line array but very narrow beam in the 
plane perpendicular to the line array) followed by processing returns 
from the entire sector, and 

b. Transmission of the acoustic pulse in a narrow beam and 
sweeping across the sector with signal processing occurring at each 
bearing increment. 

Essentially, after a suitable gate time after pulse trans- 
mission, each sampled hydrophone output provides a time series record 
of the scattered acoustic amplitudes arriving at the line array. 

Storing the results (after A/D conversion) at each sample time, n, of 
all K hydrophone outputs causes the memory to contain a time represen- 
tation of the scattered acoustic pressure field, D(n,k), k=l..K, n=l..N 
where N is the maximum sample time. The presence of a scattering point 
will be indicated by the presence in memory of signal amplitudes 
occurring at the appropriate times in each hydrophone record (a pattern 
representing the wavefront curvature). 

The sequence of time indices which uniquely represent a scat- 
tering point are available as the pre-calculated trace for that point. 

An estimator of the presence of a pattern is the summation of the signal 
amplitudes at the appropriate time indices or addresses in each of the 
hydrophone records. For example, if the signal amplitude at the 
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wavefront (from a point at r,0) was considered unity, then the total 
summation along n.(r,6) would equal K, the total number of sensors: 



As a result of noise and errors in geometry the sum would be less than 
K but would nevertheless have a maximum if a scattering point was 
present. 

In summary, what is required is to form in memory a represen- 
tation of the scattered acoustic field, search through it systematically 
for wavefront patterns by summation along pre-calculated traces for the 
ranges and bearings of interest, and recording the resultant total 
amplitudes as a function of range and bearing. A decision process 
would then follow to determine which amplitudes represented point 
scatterers, as opposed to solely noise, and then these could be dis- 
played forming the acoustical image of the insonified region in terms 
of range and bearing or rectangular coordinates. 

2. Data Structures 

To visualize the signal processing operation it is useful to 
discuss it in terms of the data structures involved: the "Indata" 
array, the "Trace" matrix, and the "Amplitude" array as shown in 
Figure 4. 

The Indata array contains the time record of the scattered 
acoustic field. It can be visualized as two dimensional in memory: 



K 




where n^ are elements of n.( r , 6 ) . 
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i.e. D(n,k) where n is the time index, and k is the hydrophone number. 
Each row would be formed by the simultaneous sampling of all hydrophones 
at intervals equal to the pulse width of 50 microseconds, beginning 
after a gate time of 6 milliseconds, and storing the results. This 
would continue for 167 samples until the last possible arrival time of 
a signal of interest (time indices 0 to 166). 

The Trace matrix, T(9,r,k) contains the traces £(r,0) for 
points in thd area of interest as identified by the intersections of all 
range and bearing increments. It can be visualized as three dimen- 
sional where the first dimension is that of bearing, the second is that 
of range, and the third is the hydrophone number. For a particular 
range and bearing, the elements in the third dimension contain the time 
index sequence n(r,9) that is the pre-calculated trace: i.e. n k = 
T(9,r,k). 

The dimensions for the Amplitude array are range and bearing. 
Each element, A(r,9), is the total sum of the contents of those 
elements from the Indata array identified by the trace £(r,9): 

K 

i.e. A(r ,9) D(n k ,k) 

k*l 

where n^ are elements of njr,9) 

Scattering points in the acoustically scanned region will be evident 
by larger amplitude elements in A(r,9) and the distribution of these 
will comprise an image of any objects present. 
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III. THE DEVELOPMENTAL ALGORITHM 



A. INTRODUCTION 

The signal processing algorithm was initially developed using a 
Hewlett-Packard 9845B desktop minicomputer (Figure 5). This machine 
uses HP enhanced BASIC [4] and was equipped with 192K of storage, two 
peripheral flexible disk drives (HP9885M and HP9885S), a Graphics ROM 
(98437B), and a CRT Graphics Memory (98470B). This proved to be a most 
convenient tool for the development of the algorithm, evaluation of 
the influence of the major parameters involved, and for the graphical 
display of simulated objects, all of which led to major conclusions and 
decisions as to the form of the algorithm and the possible physical 
acoustic imaging system. 

B. DESCRIPTION OF THE ALGORITHM 
1 . Parameters 

The parameters involved in the algorithm (Appendix A) were 
those which describe the physical implementation of an underwater 
acoustic imaging system and the area or sector of interest within the 
sediment that is to be insonified. These parameters were varied to 
evaluate their influence on the performance of the acoustical imaging 
system. 

The lengths of the horizontal line array used in the develop- 
mental algorithm were 10 meters and 5 meters; indicative of the desire 
for the largest aperture possible for angular resolution, and the most 
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Figure 5. The Hewlett-Packard 9845B 






probable realistic maximum for a physical structure to be attached to 
a maneuvering underwater vehicle. Equally spaced hydrophones along 
the line array numbering 16 and 32 were used, as were pulse widths of 
100 and 50 microseconds. The sound speed was taken as 1500 meters 
per second throughout the transmission path of the acoustic pulse. 

The sector of interest within the sediment was defined as 
lying between the ranges of 5 and 10 meters from the center of the line 
array and between bearings of -30 and +30 degrees from the normal to the 
line array. Range increments of 0.1 meters, and bearing increments of 
two degrees were chosen. This corresponded to 51 range indices from 0 
to 50, and 31 bearing indices from -15 to 15. Bearing increments of 
one degree were used in the evaluation of array bearing resolution. 

2. Calculation of Traces 

The trace, ^(r,9), for a scattering point at r,6 consists of 
the time indices representing the arrival of the scattered wavefront 
at each of the sensors in the line array. With the projector at the 
center of the line array, the distance from the projector to each 
sensor was calculated. A simple trigonometric relation was used to 
calculate the travel time for an acoustic pulse from the projector to 
a point P(r,0) and then to the k'th sensor. 

The geometry of two cases is illustrated in Figure 6 and 
Figure 7. The first case, Figure 6, applied when the k'th sensor was 
past the array midpoint; i.e. k > K/2 where K is the number of sensors 
in the line array. For this situation the trigonometrical relation 
providing the distance d^ from the k'th sensor to the point P(r,9) in 
terms of the range r, bearing 0, and sensor to projector distance b^, is 
as fol 1 ows : 
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Figure 6. The geometry for k > K/2. 





Figure 7. The geometry for k $ K/2. 
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d k 2 = t> k 2 + r 2 - 2 r cos(9O-0) 

Similarly, for the case k £ K/2 (Figure 7), the distance d k is 
obtained from 

d^ 2 = b^ 2 + r 2 - 2 b k r cos(90+8) 

The total travel distance for a scattered pulse, r + d^, divided by 
the sound speed, gives the time that the leading edge of the pulse 
arrives at the k'th sensor after transmission. The period of 6 milli- 
seconds is subtracted from the travel time, which is equivalent to 
starting data acquisition after gating out the first 4.5 meters in 
range from the center of the line array. 

The integer portion of the quotient of this time and the 
sampling interval (equivalent to the pulse width) provides the time 
index n k . Computed for each of the K sensors this results in the 
sequence n/r,8) = {n|<)- where k = 1..K. This computation was repeated 
at all ranges for each of the positive bearings. Symmetry about the 
normal to the line array permits the calculations of travel times for 
positive bearings to suffice. Negative bearings use the mirror image 
of the trace sequence; i.e. n_(r,-0) * {n^} where k = K. .1 . 

3. Simulation of Objects 

In order to develop the algorithm and in the absence of real 
hydrophone inputs, sampling, and I/O operations it was necessary to 
simulate the Indata array D(n,k). An idealized situation consisting of 
zero noise and unity signal amplitude at the wavefront of a scattered 
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pulse was considered. In other words no attempt was made to model a 
noise environment or the reflectivity of a buried object. 

A point scatterer was simulated by using the pre-calculated 
trace, as stored in the Trace matrix T(9,r,k), to identify the time 
indices n of D(n,k) that represented the wavefront arrival at each 
sensor. Unity signal amplitude was inserted into these elements. For 
the narrow acoustic beam case only a single point scatterer was assumed 
to exist at each range increment along the bearing of the beam. It was 
more complicated in the case where the entire sector was assumed insoni 
fied. In this case the acoustic field would be composed of the inter- 
ference pattern of numerous wavefronts if the sector contained a 
distributed object composed of a number of point scatterers. 

For the latter case two input situations from the sensors were 
simulated: a nonlinear input, and a linear input. The nonlinear 

situation would arise if the outputs from the hydrophones were ampli- 
fied, passed through a thresholding comparator circuit, and a binary 
input presented to the computer I/O. The unity amplitude mentioned 
above simulated this situation. If the sampled outputs from the 
sensors were subjected to A/D conversion prior to storage in memory, 
then a more accurate or linear representation of the actual pressure 
field would be available. This was simulated by the summation of the 
wavefront amplitudes in the appropriate elements of D(n,k) as each of 
the points comprising a distributed object was considered. 

The angular resolution of the imaging system was evaluated by 
simulating discrete points in the area of interest. Imaging perform- 
ance was studied using simulated objects composed of a number of 
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scattering points. The resulting Indata arrays were then operated 
upon by the algorithm to form the Amplitude array which represented the 
acoustic image of the insonified sector. 

4. Forming the Amplitude Array 

The magnitude of each element in the Amplitude array, A(r,9), 
is an estimator of the presence of a point scatterer at the range r, 
and the bearing 0. Each element results from the sum of those elements 
in the Indata array, D(n,k), identified by the trace njr,9) stored in 
the Trace matrix, T(9,r,k). 

For the case of the entire sector being insonified, the Ampli- 
tude array was formed by summing along the traces for all ranges at 
each bearing in turn for the entire sector. For the narrow beam case, 
only the traces corresponding to all ranges at the beam bearing were 
utilized, followed by the acquisition of a new Indata array for the next 
bearing. Then the traces for all ranges at that bearing were used to 
guide the summation process, continuing until the narrow beam had 
scanned the entire sector and the contents of the Amplitude array 
represented the presence of any point scatterers within the sector. 

5. Graphics Display 

To evaluate angular resolution and imaging performance of the 
proposed system the graphics capability of the HP9845B was utilized. 

It was possible to draw graphs representing the received amplitudes 
from point scatterers and thereby obtain the system's angular resolu- 
tion. Both a polar or range-bearing grid and a rectangular grid were 
devised to provide a means of referencing a displayed image. It was 
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decided that the rectangular grid would be more practical and 
meaningful to an operator of an acoustic imaging system. 

Due to the choice of 0.1 meter increments in range and two 
degree increments in bearing, a single point scatterer was represented 
by graphically "filling-in" an area on the screen equivalent to these 
dimensions and centered at the range and bearing of the point. This 
was termed a "target" in the developmental algorithm. A threshold 
was used to determine when the magnitude of an element in the Ampli- 
tude array indicated a target present. An object was simulated by 
using a number of these discrete points which the graphics display 
then presented as a number of adjacent targets. This gave the image 
outline a stepped appearance rather than the smooth lines of a real 
object. 

C. RESULTS 

1 . Range Resolution 

The use of an extremely short pulse of 50 microseconds has 
inherent range resolution much greater than required or utilized in the 
development of the algorithm. It was decided that a range resolution 
of 0.1 meters was adequate to localize and classify large buried 
objects. Accordingly, it was the range increment chosen for the pre- 
calculation of all the traces, and the subsequent processing and 
display of the acoustic image. An additional consideration was the 
memory size available to store the Trace matrix. 
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2 . Bearing Resolution 



For the case where the entire sector is insonified by one pulse, 
the effects of a number of parameter changes on the angular resolution 
of the acoustic imaging system were obtained. The bearing resolution 
was determined by the ability of the system to distinguish between two 
discrete points at the same range but at slightly different bearings. 

The two points were simulated as being in the vicinity of a bearing of 
zero degrees, and then a bearing of thirty degrees, and at ranges of 
five meters and ten meters. 

The algorithm processed the resulting Indata array in each 
case and the bearing resolution was determined from the resulting 
amplitude distribution in bearing at the range of the points. Figure 
8 gives an example of the amplitude distribution for two points at a 
number of bearing separation angles. The amplitude distributions were 
superimposed as the bearing separation between the points was increased 
from three degrees to six degrees in one degree increments. The separa- 
tion in bearing that resulted in 3db and 6db dips between the peaks in 
the amplitude distribution were recorded and are shown in Table 1. 

This was done for both the nonlinear input case and the linear 
input case, and for sixteen and thirty-two hydrophones in the line 
array. It was observed that little difference in bearing resolution 
resulted as could be expected for the idealized zero noise situation 
that was being simulated. 
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Figure 8. The amplitude distributions for two points at a number of 
bearing separation angles at 10 meters range from a 5 meter 
line array. 
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TABLE 1. Bearing Resolution as a function of 
system parameters 



3. Simulated Images 



Three large objects, each considered as being composed of a num- 
ber of point scatterers, were simulated in the sector or the field of 
view of the system. One was small and rounded as if it were a cross 
section of a long cylindrical object. The other two were long narrow 
objects; one was in a horizontal position, and the other at an inclined 
aspect angle. Figure 9 shows the positions and orientations of the 
simulated objects. 

For the case of a pulse insonifying the entire sector and 
with nonlinear inputs from the receiving hydrophones, the maximum 
amplitude for any element in the Amplitude array was equal to the number 
of hydrophones. The threshold for the display of the resulting targets 
was varied to see the spatial distribution (in range and bearing) of 
the amplitudes of the image elements. Figure 10 shows the smearing of 
the objects resulting from the lack of angular resolutions; i.e. the 
skirts of the amplitude distribution in bearing of each scattering 
point were being detected. 

Similarly Figure 11 shows the results of having linear inputs 
from the hydrophones. The largest amplitude elements were considerably 
larger than the nonlinear case, and the contours of the amplitude dis- 
tribution in the image were much steeper. It was observed that the 
size, shape and aspect angle of the objects had a substantial effect 
on the resulting amplitude distribution in the image. This would be 
beneficial for the identification or classification function of an 
acoustic imaging system particularly if it provided an operator con- 
trolled threshold for the display or if the intensity of the display 
depicted the amplitudes of the image elements. 
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non-linear system. 
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Figure 11. Imaging of the simulated objects by the 
linear system. 



For the case of a narrow two degree beam width system that 
transmits a pulse at each bearing increment, the resulting image after 
the sector was completely scanned is shown in Figure 12. It was ob- 
served that (as expected) the image exactly duplicated the simulated 
objects. 

D. CONCLUSIONS 

The developmental algorithm was used to study some of the effects 
of the system parameters on bearing resolution, and the imaging per- 
formance of possible implementations of an acoustic system that used 
the proposed signal processing algorithm. It was concluded that the 
acoustic imaging system should be a narrow beam system that scans the 
sector by transmitting a very short pulse at each bearing increment, 
provides linear inputs to memory from the hydrophone line array, and 
provides an operator controlled threshold for the image display. Al- 
though the benefits of the latter two were not evident in the simulated 
idealized narrow beam images, it is felt that they would be significant 
in a real system where noise and objects with different scattering or 
reflective properties were present. This is also the basis for 
requiring a line array as a large receiving aperture, together with the 
proposed signal processing algorithm, since objects with specular 
reflective properties would not necessarily scatter the acoustic energy 
in the direction of the projector. 
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Figure 12. Imaging of the simulated objects by the 
narrow-beam linear system. 



IV. THE NEAR REAL-TIME ALGORITHM 



A. INTRODUCTION 

An acoustic Imaging system mounted on a maneuvering underwater 
vehicle would be required to provide an operator with a near real- 
time display of any objects buried in the sediment as the vehicle 
conducts a search pattern over the ocean floor. Although the signal 
processing algorithm that was developed was basically simple computa- 
tionally, it was nevertheless very time consuming and would have been 
inadequate in an operational environment. Consequently, it was re- 
designed for increased execution efficiency and implementation on 
multiple microprocessors while maintaining the basic signal processing 
approach that has been described. Optimization of the algorithm by 
removing time consuming operations, together with incorporation of 
parallel and pipeline processes and double buffering concepts were the 
approaches taken. Experiments were conducted using a number of single 
board computers to test and demonstrate the feasibility of the rede- 
signed algorithm for near real-time execution. 

B. THE TIME PROBLEM 

An underwater vehicle acoustically scanning the sediment and moving 
at 1 meter per second would probably require a complete sector scan at 
least once per second as a minimum adequate display rate; i.e. a two 
degree beam width, approximately 0.34 meters wide at 10 meters range, 
would scan a 60 degree sector perpendicular to the direction of travel 
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once every meter of forward travel. At each bearing increment data 
acquisition would begin 6 milliseconds after pulse transmission and 
continue for 8.5 milliseconds. Thus the processing of the Indata 
array can begin 14.5 milliseconds after pulse transmission, and would 
need to be completed in about 15 milliseconds for each of the bearing 
increments if a scan were to be completed in about one second. It was 
obvious that the developmental algorithm would require considerable 
changes for an efficient microprocessor implementation. 

C. ALGORITHM OPTIMIZATION 

The algorithm was first examined with the objective of removing 
time consuming operations such as multiplication. These occurred in the 
addressing of the multi -dimensional arrays in the algorithm which were 
stored in row major order. Four additional one dimensional arrays were 
introduced which provided a table lookup of the row addressing. For 
example, the addresses for the bearing dimension of the Trace matrix, 
T(9,r,k), were stored in B_array(b), and those for the range dimension 
in R_array(r). Instead of the addressing for a single element of 
T(0,r,k) requiring three additions and two multiplications, it was 
addressed using T(B_array(b) + R_array(r) + k); i.e. five addition 
operations. This was done for each of the multi -dimensional arrays, the 
Indata array D(n,k), the Trace matrix T(0,r,k), and the Amplitude array 
A(r,9), thus reducing the data structures in the algorithm to strictly 
one dimensional arrays. 

With the intention of testing the feasibility of a microprocessor 
implementation of the algorithm using a number of Intel SBC 80/20 Single 
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Board Computers (Intel 8080 CPU) [ 5 ] , the algorithm (Appendix B) was 
written using the high level language PL/ 1-80 [6^ . Executing the 
program on one SBC 80/20 it was found that after the pre-calculation of 
the traces in T(0,r,k) the signal processing aspect of the algorithm 
required 6.45 seconds to complete one sector scan (approximately 210 
milliseconds at each bearing, and not including the time required for 
pulse propagation and data acquisition). 

D. ALGORITHM PARTITIONING 

The structure of the signal processing algorithm was examined in 
order to identify functions that could be accomplished by parallel and 
pipelined processes executing on a number of microprocessors for in- 
creased time efficiency. The concept of double buffering was also 
incorporated to reduce execution time. 

1 . Parallel Processing 

The processing carried out at each of the range increments on 
the Indata array, for the bearing of the transmitted pulse, is indepen- 
dent of any of the other ranges. This is the summation that results 
in an element of the Amplitude array A(r,0). It was decided to parti- 
tion this function into separate processes that were identical except 
for the scope of the range indices. These were called the Sum_ampli- 
tude processes and were executed in parallel on a number of micropro- 
cessors. Using three microprocessors for example, process S-| performed 
summations along the traces for the range indices 0 to 16, process S 2 
used range indices 17 to 33, and process covered range indices 34 to 
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50. In this manner all ranges were processed in one third the time 
that was required by a single process executing all range increments. 

2. Pipeline Processing 

The function of data acquisition is called the Indata process 
and was seen to be amenable to pipeline processing with the Sum_ampli- 
tude processes. This process could control the timing of pulse trans- 
mission, sampling the hydrophone outputs, and storing the results into 
two buffer arrays. The Indata array of the developmental algorithm 
became two arrays called the A_indata array and the B_indata array. 
After the A_indata array was filled with the data from a particular 
bearing, the Sum_ampl itude processes were permitted to begin their 
operations on the A_indata array. Simultaneously, the Indata process 
proceeded to fill the EMndata array with the data from the next 
bearing increment. When the Sum_ampl itude processes had finished with 
the A_indata array, they processed the EMndata array. The Indata 
process could then utilize the A_indata array for the next bearing and 
so on in a pipeline fashion. Figure 13 illustrates the relationship 
between the processes. 

3. Semaphores 

Communication between the processes was accomplished with 
binary semaphores. The Indata process set the semaphores A_data and 
B_data whenever it was utilizing either the A_indata or B_indata array. 
Similarly, each of the parallel Sum_ampl itude processes set a semaphore 
depending on the data array it was processing. Thus for three parallel 
Sum_ampl itude processes six semaphores were utilized: A_ampl , A_amp2, 
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gure 13. The relationship between the processes. 



A_amp3, and B_ampl , B_amp2, B_amp3. The Indata process waited for each 
of the Sum_ampl itude semaphores for one of the data buffer arrays to be 
unlocked before it could. proceed to fill the buffer with new sampled 
data from the hydrophone outputs. Similarly, each of the Sum_ampli- 
tude processes waited on the Ajdata or B_data semaphore before it 
began processing the data. In this way completely independent but 
synchronized operation of all the processes was achieved. The opera- 
tional form of the algorithm in Appendix C illustrates the parti- 
tioning of the processes and the use of the semaphores. 

E. EXPERIMENTAL DEMONSTRATION 

1 . Equipment 

To test and demonstrate the feasibility of near real-time 
execution of the redesigned algorithm an Intel lec 800 Microcomputer 
Development System (MDS) [7] , shown in Figures 14 and 15, with the 
CP/M monitor control program [8] was used. The MDS is capable of 
being used as a partial system simulator to check out software execu- 
ting from RAM in the Intellec MDS. A total of 64 kilobytes of memory 
was available and four Intel SBC 80/20 Single Board Computers [5] were 
inserted into the MDS motherboard as shown in Figures 16 and 17. 

2. The Processes 

The contents of the Trace matrix, the pre-calculated traces, 
were computed using the program of Appendix B, and stored in a floppy 
disk data file for subsequent loading into memory. The redesigned 
"operational" algorithm (Appendix C) was written as if it were to be 
executed on a single processor, and consisted essentially of a Setup 



41 




42 



Figure 14. The system used for the experimental demonstration 
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Figure 15. The Intel lec 800 Microcomputer Development System (MDS) 







WMM 
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Figure 16. The Intel Single Board Computer SBC 80/20 
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inserted in motherboard. 





procedure that calculated the contents of the address lookup arrays, the 
Indata procedure, and the Sum_amplitude procedure. This was compiled and 
linked so the program started at location 3000 Hex. The choice of 3000 H 
as the initial address for the program was based on the fact that addres- 
ses 3000 H to 3EFF H were available as private RAM on each of the SBCs. 
All the data structures and variables such as the semaphores that were 
to be global to all the processes were located in common memory above 
4000 H that was accessible by all the SBCs. An empty dummy array of the 
necessary size was declared to ensure the assignment of common memory 
addresses to the global variables. 

Using the MDS and the Dynamic Debugging Tool (DDT) program [9] 
the program was altered at the machine code level (hex representation) 
by removing the appropriate call statements and creating two versions: 
one which executed only the Indata procedure, and the other the Sum_am- 
plitude procedure. In this way the two processes stored in the private 
memory of separate SBCs could operate independently but in a syn- 
chronized manner using the semaphores in common memory. The data that 
was operated on, as well as all the data structures, were available in 
common memory. 

The Sum_ampl itude version was then amended using DDT so that all 
the local variables in the procedure, which the compiler had assigned to 
memory locations in the common area, were readdressed so they were 
within an SBCs private memory. The locations assigned to the dummy 
array were utilized for this purpose. This then provided a template for 
three subsequent copies of the Sum_ampl itude process each of which was 
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individually altered to process a segment of the range indices, and then 
stored on a disk file. 

Similarly, the addresses of variables that were local to the 
Indata procedure were amended so they were within an SBC's private 
memory. Having loaded the Trace matrix from the floppy disk data file 
into the appropriate memory locations that had been assigned by the 
compiler, the entire Indata version (which included the data structures) 
was saved on a disk file. A subsequent reloading of this file into 
memory caused the Trace matrix to also be available and so a separate 
loading of the Trace matrix data file was not necessary. For the 
experimental demonstration of the feasibility of the multiple micropro- 
cessor concept the Indata process controlled the bearing parameter and 
caused the entire algorithm to scan ten complete sectors and then stop. 
An operational version would continue indefinitely. As the double buf- 
fering design resulted in two bearing increments being processed on each 
pass through the procedure, it was necessary to ensure that each sector 
consisted of an even number of bearing increments, i.e. thirty vice 
thirty-one. No attempt was made to include pulse propagation time or 
data acquisition time. The procedure that would fulfill these functions 
was left as a "stub" since it was not required for the demonstration of 
the synchronization between the processes. 

3. Results 

The capabilities of the MDS and DDT were used to subsequently 
load the Indata process into the private memory of one of the SBC's, and 
the Sum_ampl itude processes into the private memory of each of the three 
remaining SBCs. Figure 18 provides a block diagram of the experimental 
configuration where process D is the Indata process and processes S-j to 
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Figure 18. Block diagram of experimental 
configuration. 
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S 3 are the Sum_ampl itude processes. With the processes executing simul- 
taneously the time required for the amplitude summation function to be 
performed on ten sector scans was recorded. 

Executing the algorithm with a different number of Sum_ampli- 
tude processes in parallel, while ensuring that all range indices were 
processed, permitted an evaluation of the effect of increasing the number 
of parallel processors. The resulting times for one sector scan and 
consequently the performance in terms of scans per second are given in 
Table 2 and Figure 19. 



Processes 


Time for one scan 
(seconds ±0.01 ) 


Scans per 
Second 


D, S 1 


6.45 


0.155 


D, S 1§ S 2 


3.23 


0.310 


o, s 1§ s 2 , s 3 


2.15 


0.465 


TABLE 2. Results of the experimental demonstration of the 
feasibility of implementing the signal processing 



algorithm on multiple microprocessors. 

The correct operation of the algorithm together with the per- 
formance times indicated that the design was a feasible approach to 
obtaining near real-time execution of the proposed signal processing 
algorithm. The fact that the processes were stored in the SBCs private 
memory and required access to the bus only when the common data struc- 
tures and semaphores were referenced resulted in a minimum of bus usage. 
An examination of the code for the innermost loop of the Sum_amplitude 
process (the summation of the Indata array elements identified by the 
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NUMBER OF PARALLEL PROCESSORS 



trace) indicated that the bus usage was 7.9% of the loop's execution 
time. The number of parallel SBCs executing simultaneously could theo- 
retically be increased to twelve before bus contention would adversely 
influence the overall execution time. Table 2 and Figure 19 depicting 
the results of the experiment clearly indicate the linear relationship 
between performance, in terms of scans per second, and the increasing 
number of parallel processors used. If the number of SBCs operating in 
parallel were increased to seven, plus one for the Indata process for 
a total of eight, then a scan rate of one sector per second would be 
achieved . 
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V. CONCLUSIONS 



For the problem of locating and classifying objects buried in 
marine sediments an underwater acoustic imaging system has been pro- 
posed and the thrust of this thesis has been the design of the signal 
processing algorithm that would extract the information from scattered 
acoustic signals for subsequent operator display. Using a developmen- 
tal form of the algorithm it was determined that the acoustic imaging 
system should be a narrow beam system for optimum resolution and 
imaging performance. Since the processing time was excessive, the 
algorithm was redesigned for implementation on multiple processors for 
execution efficiency. It was experimentally demonstrated that this 
design was a feasible approach to obtaining near real-time execution of 
the algorithm. 

The maximum rate of sector scanning that is possible for the pro- 
posed system is slightly in excess of two sector scans per second. 

This was arrived at by considering the pulse propagation time and the 
data acquisition time, a total of 14.5 milliseconds at each bearing, 
as the minimum time between pulses. It has been shown that the use of 
multiple microprocessors can result in an efficient implementation of 
the proposed signal processing algorithm. The number of parallel pro- 
cessors required to reduce the amplitude summation process to about 15 
milliseconds (and thus achieve the maximum rate of two scans per 
second) would be dependent on the microprocessor chosen. 
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VI. RECOMMENDATIONS 



The signal processing algorithm and its implementation on mul- 
tiple microprocessors, as discussed in this thesis, would form the core 
of a special purpose processor for the acoustic imaging system. Two 
areas would require development and detailed definition before the 
entire special purpose processor could be designed: data acquisition 

to provide inputs to the processing algorithm, and output from the 
algorithm in the form of a display. 

It is recommended that the Indata process outlined in the thesis 
be expanded from the current “stub" status to perform the functions 
previously mentioned; i.e. control of pulse transmission timing, 
sampling the hydrophone outputs, and storing the results into the two 
buffer arrays. It is anticipated that the approaches used for the 
algorithm design would be readily applicable to these functions. How- 
ever, the actual details would be highly hardware dependent which would 
require investigation prior to designing a parallel/pipeline approach to 
this portion of the specialized processor. 

Similarly, the interface between the multiple microprocessor 
system (performing the signal processing) and a display processor would 
require definition and would be hardware dependent. It is considered 
likely that the image information after each sector scan could be 
stored thus forming a three dimensional representation in memory (the 
third dimension resulting from the travel of the underwater vehicle 
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over a section of sediment). This would provide an opportunity for a 
two dimensional display in either the vertical plane (such as range and 
bearing) or in a horizontal plane at some range. Furthermore, the possi- 
bility of a three dimensional display could be considered. It is evi- 
dent that there is much promising work to be accomplished in order to 
produce the special purpose processor for the acoustic imaging system. 
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3190 1 PLOT 0,-<Max_rar.cie + l> 
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3560 FOR Y=-Max_range TO -M i rt_rar>ge + 1 
3570 MOVE -4 . 9 , Y 

3580 LABEL Y 

3590 NEXT Y 
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3530 FOR Bear i ng=- 1 Q TO -8 STEP 2 
3540 FOR Range=5 . 2 TO 5.4 STEP 

3550 GOSUB Set up_ ob j ec t s 

3560 NEXT Range 

3570 NEXT Bearing 



3980 FOR Bear i ng = -6 TO -4 STEP 2 
3990 FOR Range =9. 1 TO 9.3 STEP 

40O0 GOSUB Setup objects 

4010 NEXT Range 

4020 NEXT Bearing 
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4339 FOR Range=7. 1 TO 7.3 STEP 

4340 GOSUB Setup_objects 

4350 NEXT Range 
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APPENDIX B. THE ALGORITHM IN PL/I -80 
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; FILL AMP_ARRAY WITH ADDRESSES CORRESPONDING 
TO THE BEARING DIMENSION OF THE AMPLITUDE MATRIX V 
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